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Abstract 

The Stern-Gerlach experiment is the fundamental experiment in order to exhibit the quantization 
of spin and understand the measurement problem in quantum mechanics. However, although the 
Stern-Gerlach experiment plays an essential role in the teaching of quantum mechanics, no complete 
analysis of this experiment using Pauli spinors is presented in the pedagogical literature. This paper 
presents such an analysis and develops implications for the theory of quantum measurement. 

We first propose an analytic expression of both the wave function and the probability density in 
the Stern-Gerlach experiment. Our explicit solution is obtained via a complete integration of the 
Pauli equation over time and space. The probability density evolution describes a slipping of the 
wave packet into two separate packets due to the measurement device, but it cannot account for 
impacts. 

We therefore calculate the de Broglie-Bohm trajectories, which not only explain impacts naturally, 
but also accounts for the spin quantization following the magnetic field gradient. It is then possible 
to propose a clear explanation of measurement effects in the Stern-Gerlach experiment. 
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I. INTRODUCTION 



As they were studying the deviation of a silver atoms beam in a highly inhomogeneous 
magnetic field (cf. FIG. [TJ Stern and Gerlach (1922)^ found empirical results which chal- 
lenged common sense prediction. Instead of being scattered, the beam split into two symetric 
beams, which produced two spots of equal intensity on a screen, at equal distances from the 
axis of the original beam. 

z 
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FIG. 1: Schematic configuration of the Stern-Gerlach experiment. 

This experiment motivated the introduction of spin quantization as an intrinsic magnetic 
moment. It also clearly exhibits the measurement problem in quantum mechanics, which 
remains an active field of study.— 

This paper brings new elements, such as the analytic expression of the wave function and 
the probability density in the Stern-Gerlach experiment. The explicit solution is obtained 
via a complete integration of the Pauli equation over time and space. As far as we know, the 
analytic presentation of the Stern-Gerlach experiment in text-books^ is only given in time, 
not in space. The first explicit calculation in space and time of Stern-Gerlach experiment 
was given by Dewdney et a^, but it remains incomplete, as is also incomplete their explicit 
solution in space and time of the Dirac equation.— Recent presentations in space and time 
are only given with numerical simulations.-^ 

The analytic solution presented here gives an explicit expression of the probability den- 
sity's evolution in space, explaining the semi-classical presentation and showing the wave 
function separation. However, this continuous evolution in space of the wave packet into 
two wave packets does not account for the particle impacts. We therefore also calculate 
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the de Broglie-Bohm trajectoriesS as we formerly did in the case of Young's double slit 
experiment.— These trajectories not only provide a natural explanation for the impact of 
particles, but also describe the spin quantization along the z-axis. It is then possible to 
propose a clear explanation of measurement in the Stern-Gerlach experiment. 

The explicit solution in terms of Pauli spinors and the evolution of the probability density 
for the Stern-Gerlach experiment are presented in section 2. The de Broglie-Bohm trajec- 
tories, as defined by the interpretation of the Pauli equation by Bohmfii and Takabayasi,— 
are simulated in section 3. Details of the calculations are provided in Appendix A. 

II. THE PROBABILITY DENSITY CALCULATION IN THE STERN-GERLACH 
EXPERIMENT 

Silver atoms contained in the oven E (Fig.HJ are heated to a high temperature and escape 
through a narrow opening. A second aperture, T, selects those atoms whose velocity, vq, 
is parallel to the y-axis. The atomic beam crosses the gap of the electromagnet Ai before 
condensing on the screen. Pi . The magnetic moments of the silver atoms before crossing the 
electromagnet are oriented randomly (isotropically). In the beam, we represent the atoms 
by their wave function ; one can suppose that at the entrance to the electromagnet, Ai, and 
at the initial time t = 0, each atom prepared can be described by a Gaussian spinor in x 
and z: 



The variable y will be treated in a classical way with y = vt. For a silver atom, one has 
m = 1.8 X 10-25 kg, vo = 500 m/s(with T=1000°K), ao=10-^m. 

In (QJ), 6q and (po are the polar angles characterizing the initial orientation of the mag- 
netic moment, 9o corresponds to the angle with the z-axis. This initial orientation being 
randomized, one can suppose that 9o is drawn in a uniform way from [0, vr] and that ipo is 
drawn in a uniform way from [0, 27r]. In this way, we give a model of a mixture of pure 
states. 

The evolution of the spinor ^ = (^■'') in a magnetic field B is then given by the Pauli 




(1) 
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equatio] 
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"Mh-.^bJ*-! (2) 

J'- J \'^- 

where fj^B = is the Bohr magneton and where cr = (a^^, dj^, a^) corresponds to the three 
Pauli matrices. 

The particle first enters an electromagnetic field B directed along the z-axis, 5^ = B'^x, 
By = 0^ B, = Bo - B'qZ, with Bq = 5 Tesla and B'q = \^\ = 10^ Tesla/m over a length 
Al = 1 cm. Such a vector B satisfies Maxwell's equations, since dif (B) = 0. 

Once it leaves the magnetic field, the particle travels freely until it reaches the screen. 
Pi, placed a distance D = 20 cm beyond the magnet. 



A. The Probability Density in the Magnetic Field 



The particle stays within the magnetic field for a time At = -^ = 2xl0 ^s. During this 
time [0, At], the spinor is (see Appendix A) 











yR.e'^ J 





cos|i(27ra2)-ie' 



' 2m ' ' 



iL,3_ 



(3) 

Since the initial spinor direction is random, the atomic density, p(z, t) is found by in- 
tegrating i?^ + R^_ on (6*0, 9?o) and x (notice that R\ + R^_ is independent of ipo). So one 
gets: 



Ms£q,2)2 



(,+ MS£at2)2 



"0 + e 



(4) 



B. The Probability Density after the Magnetic Field 



After the magnetic field, at time t + At (t > 0), the spinor becomes (see Appendix A) 



where 



u = = Im/s. 



2m m 
One can deduce (as previously) the atom density p at (z,t + At): 



p{z,t + At) = (27ro-, 
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(z — —ut)'' 
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(6) 



(7) 



Figure [2l shows the probabihty density of the silver atoms as a function of z at several 
values of t ( the plots are labelled with y = vt). The beam separation does not appear at 
the end of the magnetic field (1 cm), but 10 cm further along. 
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FIG. 2: Evolution of the probability density of silver atoms. 



In equations (|H) and (|Z|) one recognizes classical trajectories: 

2^ = |±/f£^ for t G [0,At];±ZA±Mt for At + t(t>0)} (8) 
m 

corresponding to particles with magnetic moments of pz = Pb and pz = —pB respectively. 
The trajectories are parabolic inside the magnetic field and linear after. The two spots 
and N~ become separated when — Z- ^ 4(To, which occurs at the separation time 

Statistically, everything happens as if the atom's moment were quantified in two parts: 
one half of the particles having Pz = I^b and the other half having Pz = —I^b- This re- 
sult explains why the semi-classical description of this Stern-Gerlach experiment is usu- 
ally used. This semi-classical description starts from the quantization of spin and de- 
duces from Ehrenfest's Theorem the average trajectories of the spinors with initial spinors 
= (27r(T^)-2e ^"0 (J) and = (27r(j^)-2e ^"o (°) respectively. Note that this is 
a statistical interpretation, and the individuality of the atoms represented by the angles 
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^0 and (fio is lost. However experimentally, one does not observe directly the wave func- 
tion, but individual impacts of the silver atoms on the Pi screen. The usual explanation of 
these individual impacts on the screen is decoherence,-^"^ caused by the interaction with the 
measurement device. Here, the evolution of the probability density given by equations (jH) 
and (IH) correctly describes a separation of the wave packet into two packets thanks to the 
measurement device, but cannot describe the individual positions of these impacts. 

III. IMPACTS AND QUANTIZATION EXPLAINED BY TRAJECTORIES 

To explain individual impacts, we simulate the silver atom trajectories in the de Broglie- 
Bohm interpretation just as we didl^ for the neon atoms in the Young's double slit ex- 
periment. In this unusual presentation of the Quantum Mechanics results, the particle is 
represented not only by its wave function, but also by the position of its center-of-mass. 

Indeed at the first instant, the wave function z) gives the initial probability density 

z). This density, which doesn't depend on h, is the classic presence density of a silver 
atom. So in classic mechanics, one does have an undetermination on the atom position, and 
in order to describe its evolution, it is necessary to precise its initial position. The principle 
of the De Broglie-Bohm interpretation is to do the same in quantum mechanics. 

So, in the Pauli equation case atoms have trajectories that are defined by using the 
center-of-mass velocity v(a;, y, z, t) give n^^'^^ by : 

h h 

2mp ^ ' 2mp ^ ^ ^ ^ 

where \I'"I' = [ip^yipl]. 

Let us show how this interpretation gives the same results as the Copenhagen school. 
One verifiesiiii^ that with v given by (fTn|l . the probability density p{x,y,z,t) = \[^^\[^ = 
I/, 2;, t)p of the \l/ spinor solution of the Pauli equation Q satisfies the Madelung 
continuity equation: 

dp 

— + div[p\) = 0. 

One can deduce from it, that if a particle family with the initial probability density po{x, y, z) 
follows the de Broglie-Bohm trajectories, its probability density at the t time, will be 
p{x,y,z,t). 
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Thus, these two interpretations give statistically identical results, but the de Broglie- 
Bohm interpretation predicts the position of individual impacts as well. We shall see that 
these trajectories also explain the spin quantization following the magnetic field gradient. 

In equation (fTT1|l . the last term ^^roti^^"^ a'if) corresponds to the Gordon current. Its 
contribution to velocity is here negligible. We will therefore not take it into account from 
now. 

In the de Broglie-Bohm interpretation, the individual particle is not only described by 
its wave function, but by its initial position (xq, zq, yo = 0) as well. So, trajectories in x and 
z are given by the differential equations: 

dx _ 1 d{S+ + S^) 1 d{S+-S-) ^^^^ 
dt 2m dx 2m dx 

dz 1 diS+ + S^) 1 d{S+-S\ , , 

— = - + ^cos^ (12) 

dt 2m oz 2m oz 

with tan | = 

A silver atom with a polarization (^o,V^o) and a position zq at the entrance of the electro- 
magnet Ai will satisfy the differential equation in the [0, At] period: 

dz fiBB'nt ^, , . , e(z,t) 00 -^^f^ , , 

— = cose(z,t) with tan ^ ' ^ = tan — e "'""o (13) 
dt m ^ ' ' 2 2 ^ ^ 

and for the At + t (t >0) period: 

dz tanh(^^^^) + cos^o ^ e(z,t) Oo -if^^ 

— = u r — r-r\ cind tan = tan — e ''o . (14) 

dt i + tanh(^^^^)cos0o 2 2 ^ ^ 

Figure El presents a plot in the x,y plane of a set of 10 silver atom trajectories with the 
initial polarization (6'o = f , V^o = 0) and whose initial position zq has been randomly chosen 
from a Gaussian distribution with standard deviation uo. The spin orientations 9{z,t) are 
represented by arrows. 

Figure 01 presents, for a silver atom with the initial polarization (6*0 = f , v^o = 0), a plot 
in the x,y plane of a set of 10 trajectories whose initial position Zq has been randomly chosen 
from a Gaussian distribution with standard deviation ctq. The spin orientations 9{z,t) are 
represented by arrows. 
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FIG. 3: Ten silver atom trajectories with initial polarization {9q = |) and initial position zq; 
Arrows represent the spin orientation 6{z,t) along the trajectories. 

One can notice that the final orientation, obtained after the separation time ts, will 
depend on the initial particle position zq in the wave packet and on the initial angle 6*0 of 
the atom magnetic moment with the axis z. 

We obtain +f if ;zo > and -f if < with 

/°=aoF-i(sin2|) (15) 

where F is the cumulative distribution function of the standard normal distribution. 

So besides explaining the position of impacts, this simulation shows that it is possible to 
give a simple interpretation of quantization on z-axis. 

Figure IH presents a plot in the x,y plane of a set of 10 silver atom trajectories whose 
initial characteristics {9o,ipo,ZQ) have been randomly chosen: Oq and ipo from an uniform 
distribution and zq from a Gaussian distribution with standard deviation (Tq. The spin 
orientations 0{z,t) are represented by arrows. 
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FIG. 4: 10 silver atom trajectories after the electro-magnet; Arrows represent the spin orientation 
9{z, t) along the trajectories. 

IV. CONCLUSION 

It is now possible to propose the following interpretation of measurement in quantum 
mechanics. There certainly exists an interaction with the measuring apparatus as it is 
classically explained, and there exists a minimum delay necessary for the measurement 
ts — Yet we believe that this measurement and this delay do not have the meanings 
that are classically attributed to them. In the present case, the measuring apparatus itself 
gives the orientation of the spin either in the direction of field, or in the direction opposite to 
the field, depending on the position of the particle in the wave packet. The measuring delay 
is then the time which is necessary for the particle to point its spin to the final direction. 

Let us notice that in this numerical study of the Stern-Gerlach experiment we didn't 
need any of the classical postulates of measurement in quantum mechanics, such as eigen- 
values of the Hamiltonian or wave packet reduction. These two postulates may account for 
experimental results, but they do not give any idea of the transitions that lead to these 
results. Instead we only used the quantum mechanics equations (Pauli equation) and made 
only one hypothesis : the centers of mass of the atoms in the atomic beam are spacially 
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distributed according to the density given by the wave function, and follow paths that are 
compatible with the continuity equation (De Broglie-Bohm hypothesis). From this one and 
only hypothesis, we provided altogether: 

- a simple explanation of the position of the particle impacts; 

- a simple explanation of the spin quantization along the measurement axis ; 

- an simple explanation of the transition towards the Hamiltonian eigenvalues. 



APPENDIX A: CALCULATING THE SPINOR OF THE STERN-GERLACH EX- 
PERIMENT 

In the magnetic field B = (5^, 0, i?^), the Pauli equation Q gives coupled Schrodinger 
equations for each spinor component 

ih——{x, z, t) = -—V'^'ip±{x, z, t) ± fiB{Bo - BQz)iij±{x, z, t) =f ifiBB'^xip^ix, z, t). (Al) 
at 2m 

If one effects the transformatior>i^ 

tjj±{x,z,t) = exp{± ^^^^°^ )^^{x, z, t) 

equation (jAljl becomes 

ih-^{x,z,t) = -—'^'^■ip±{x,z,t) T f^BB'QZip^{x,z,t) Tif^sB'^xip^ix, z,t) exp{±i — ) 

The coupling term oscillates rapidly with frequency on = = 1,4 x lO^^s^^. Since 

I-BqI ^ \Bqz\ and |i?o| ^ l-^o^L the period of oscillation is short compared to the motion of 
the packet along its trajectory. Avering over a period that is long compared to the period 
of oscillation, causes the coupling term to vanish, yieldingi^ 

ih-^{x, z, t) = -^^^^±(2^5 z, t) =F HbBoZiP^Ix, z, t). 



The initial wave function ipj.{x, z) = ip±{x, z) = ipx{x)'^±{z) with ip'^^x) = {271(7^) ^"o ^ 
i'+iz) = {2naQ)^ie *'^o 005^6*^ and i)'^{z) = (27rcro)~2e ■*''o i sin ^e"*" allows a separa- 
tion of variables x and z. Then we can compute explicitly the preceding equations for all t 
in [0, At] with At = f = 2 x lO^s. 
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We obtain: 'iljj.{x,z,t) = ilj^{x,t)ipj.{z,t) with 



i h , 
exp-[--tan ( 



ht 



+ 



(A2) 



^+{z,t) = ipK{z,t) cos ^e^^ and K = —{jleBq 



tjj_{z,t) = ^ljKiz,t)ism * 2" and K = +^bB'q 



- — tan-V ^ - Ktz + ^ ^HiZ.^^ 1 



^K{z,t) = {271 (j^)-^e exp^[--tan-\^^) - Ktz - ^^^ + 

(A3) 



Equations ()A2j) and ()A3|) are classical results.— 

hA.t 

2mcr() 



The experimental conditions give = 4 x 10 m ^ ctq = 10 m. We deduce the 



.2 



approximations at — ctq, ^lj^{x,t) ~ ip^ix) = (27rcrQ) 4e ^ and 



~ (27ra2)^4e ^-cS exp -[-Kt^; - -^]. (A4) 
At the end of the magnetic field, at time At, the spinor equals to 

^{x,z,At)=Mx,At)\ ^^^'^'^^^ I (A5) 

^.{z,At) 



with 



At) = (27ra2)-3e cos ^e^^+ 



ij_{z, At) = {271 a^yie ^^'t ""^"Isin— e 



1 -(£±£a) irauz 



za = , u = and 

2m m 



ipo I^BBoAt K^{Atf ifo , HoBoAt K^{At) 



3 



= t; Z n ; 'Z'- = — K- + 



2 h 6mh 2 h 6mh 

We remark that the crossing through the magnetic field gives the equivalent of a velocity 
+u in the direction Oz to the function ip^ and a velocity —u to the function ip^. Then 
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we have a free particle with the initial wave function ()A5jl . The Pauli equation resolution 
gives again ip±{x,z,t) = ipx{x,t)ip±{z,t) and with the experimental conditions we have 

ip^{x,t) ~ (27r(To)"4e ^ and 

+ At) ~ (27r(T2)-4e ' cos — 
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